2013 15:36 



World Scientific Review Volume - 9in x 6in 



Chapter 1 

Superfluid Pairing in Neutrons and Cold Atoms 

J. Carlson, S. Gandolfi 

Theoretical Division, Los Alam,os National Laboratory, Los Alamos, New 

Mexico, 87545, USA 

Alexandres Gezerlis 

Department of Physics, University of Washington, Seattle, Washington, 

98195, USA 

Ultracold atomic gases and low-density neutron matter are unique in 
that they exhibit pairing gaps comparable to the Fermi energy which 
in this sense are the largest in the laboratory and in nature, respec- 
tively. This strong pairing regime, or the crossover between BCS and 
BEC regimes, requires non-perturbative treatments. We describe Quan- 
tum Monte Carlo results useful to understand the properties of these 
systems, including infinite homogeneous matter and trapped inhomoge- 
neous gases. 

1. Prom cold atoms to neutron stcirs and back 

Although the energy and momentum scales of cold atomic gases and atomic 
nuclei differ by many orders of magnitude, we can gain insight into super- 
fluid pairing in the strongly correlated regime by comparing and contrasting 
the two systems. Dilute neutron matter and ultracold atomic gases near in- 
finite scattering length (unitarity) are similar in that their equation of state 
and pairing gaps, when measured in terms of the Fermi energy, are compa- 
rable. For this reason these systems can be viewed as "high-temperature 
superfluids", even though one occurs at very small and the other at very 
large temperature. Since the pairing gap A is roughly proportional to the 
critical temperature Tc, both these systems exhibit very strong pairing, 
the strongest ever observed before in nature or experimented with in the 
laboratory. In cold fermionic atoms the particle-particle interactions can 
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be tuned experimentally, thus also mimicking the setting of low-density 
neutron matter, which is beyond direct experimental reach. More specif- 
ically, at low densities s-wave scattering can be described using only two 
parameters: the scattering length a and the effective range r^,. 

1.1. Ultracold Atomic Gases 

For ultracold fermionic atoms at low temperature,'^' the superfluid phase 
arises across the entire spectrum of weak to strong attractive interactions. 
Experiments with ^Li or "^^K have an interparticle spacing that is signif- 
icantly larger than the characteristic length of the interaction. Over the 
last decade experimentalists managed to vary the interaction strength (the 
scattering length a) across a resonance, through a regime known as "uni- 
tarity". The cold atom gases are very dilute, so the typical scale of the 
interaction range (or the effective range r^) is much smaller than the av- 
erage interparticle spacing. For these broad Feshbach resonances we have 
kpTe ^ 1, where kp is the Fermi wave vector {p = fc|,/(37r^)), so at fixed 
density the effective-range can be taken to be very small, essentially zero. It 
may be possible to use narrow and wide resonances in cold atoms to study 
the case of varying rg experiment ally,'^! and thus directly simulate neutron 
matter. 

A large variety of equilibrium and dynamic properties have been mea- 
sured in cold Fermi atom experiments. Experiments using ^Li at Duke Uni- 
versity'^ and at EN^have measured the ground-state energy of the system, 
essentially finding it to be in good qualitative agreement with Quantum 
Monte Carlo (QMC) predictionJ^'21 The ground-state energy per particle 
is conventionally given in units of the energy of a free Fermi gas at the 
same density Epc = ?>Ep/b — 3tfkp/{10m) a,s E = £,EpQ. Recent sign- 
free auxiliary field QMC calculation^[i2] and experiments^ give very precise 
results ^ at unitary, ^ — 0.372(5) and 0.375(5), respectively. 

The great advantage of cold atom systems is that they offer quantitative 
experimental results in the strong-coupling regime with simple interactions. 
An important example is the pairing gap at unitarity. Experiments at 
MIT and Rice probed lithium gases with population imbalance (also called 
"polarized" gases). An MIT experimenlJi^l established the phase diagram 
of a polarized gas, revealing spatial discontinuities in the spin polarization. 
This experiment was then usecF^ to extract the pairing gap, which was 
found to be approximately half of the Fermi energy Ep, in good agreement 
with QMC calculations:^ The gap is conventionally given in units of Ep as 
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A — TjEp. The MIT group later used RF spectroscopy to independently 
determine the gap, finding it to be in agreement with the afore-mentioned 
calculation and extractioniSl 

1.2. Neutron matter and neutron drops 

As already mentioned, at very low densities, neutron matter is very similar 
to cold Fermi atoms. The neutron- neutron scattering length is fixed and 
large (a « —18 fm), but by varying the density we can probe different 
values of kpa. Neutrons in the inner crust of a neutron star are expected to 
pair in the -'^Sq channei P^ I ^^ I at higher densities, also important in neutron 
stars, the effective range becomes important as do the repulsive parts of 
the neutron-neutron interaction, and higher partial waves. 

The inner crust of a neutron star contains a neutron gas embedded in a 
sea of ions. The properties are largely determined by the EOS and pairing of 
homogeneous matter. Pairing in neutron matter has been studied for many 
decades, leading to a large spread of predictions of the ^Sq pairing gap even 
for this idealized system. The gradient terms in the density functional and 
the behavior of the pairing gap in an inhomogeneous system are potentially 
important for the crust of a neutron star and can be studied simulating 
neutrons in external fields.'^ 

It may be possible to access superfluidity in neutron star matter observa- 
tionally: superfluidity in a neutron star is often used to explain its dynam- 
ical and thermal evolution, impacting the specific heat, bremsstrahlung, 
and pair breaking/ format ion EHHI Additionally, a cooling mechanism that 
makes use of superfluid phonon^MUl has been proposed. Whether this 
mechanism is competitive to the heat conduction by electrons in magne- 
tized neutron stars or not is a question that is directly correlated to the 
size of the gap. 

Neutron matter and neutron drop computations also hold significance in 
the context of traditional nuclear physics: equation of state results at densi- 
ties close to the nuclear saturation density have been used for some time to 
constrain Skyrme and other density functional approaches to heavy nuclei, 
while the density-dependence of the ^Sq gap in low-density neutron matter 
has recently also been used to constrain Skyrme-Hartree-Fock-Bogoliubov 
treatments in their description of neutron-rich nuclei.t^ Recent ab initio 
results for neutron drops point to a need for more repulsive gradient terms 
in inhomogeneous neutron matter, and a reduced isovector spin-orbit and 
pairing strength compared to standard functionals.l^ 
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2. BCS and Quantum Monte Carlo methods 

BCS theory has been critical to understanding many of the pairing proper- 
ties of nuclei. Quantum Monte Carlo many-body simulations, on the other 
hand, have been used for some time to calculate the equation of state of 
strongly-correlated systems, e.g. liquid helium. In such systems, however, 
QMC methods were unable to reliably calculate pairing gaps because of the 
vast difference in scale between the energy of the entire system (of the order 
of eV's per particle) and the pairing gap (of the order of meV). Thus, the 
same feature of strongly paired fermionic systems (namely the large pairing 
gap) that precludes the application of mean-field theories is precisely the 
reason that allows many-body simulation techniques to be used. 

2.1. Weak coupling 

We first briefly review the weak-coupling regime where exact results are 
available. At extremely low densities (|A:Fa| << 1) the effective coupling 
between two fermions is weak and matter properties can be calculated an- 
alytically. The ground-state energy of normal (i.e. non-superfluid) matter 
in this regime was calculated by Lee and Yang in 1957P^ 

= 1 + fc^a + (11 _ 2 In 2) {kpaf , (1) 

where Epc is the energy of a free Fermi gas at the same density as the 
interacting gas. While this expression ignores the contributions of super- 
fluidity, these are exponentially small in [l/kpo). In the next section we 
compare these results to QMC calculations for \kpa\ > 1. 

The mean-field BCS approach celebrated in the present volume reduces 
in the weak-coupling limit to: 

As was shown in 1961 by Gorkov and Melik-Barkhudarov,'^^ the BCS result 
acquires a finite polarization correction even at weak coupling, yielding a 
reduced pairing gap: 

.n., ^ 1 8 H^kl / TT \ 

^(^-)-(4^^;^^^^H^) ■ 

Thus, the polarization corrections reduce the mean-field BCS result by a 
factor of l/(4e)^/'^ ~ 0.45. Interestingly, if one treats the polarization 
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effects at the level of sophistication used in the work of Gorkov and Melik- 
Barkhudarov, this factor changes with fcj?a,^^' though there is no a priori 
reason to expect such an approach to be valid at stronger coupling [kpa 
of order 1 or more). Calculating the pairing gap in this region has been a 
difficult task, as can be seen from the multitude of publications devoted to 
this subject over the past twenty years EHMHIMI 

2.2. Quantum Monte Carlo 

Quantum Monte Carlo simulations typically begin with a local Hamiltonian 
of the form: 



where N is the total number of particles. In the case of cold atoms the 
interaction only acts between opposite spin pairs and is effectively a contact 
interaction (often simulated on the computer using a short-range potential) . 
The neutron-neutron interaction is more complicated, containing one-pion 
exchange at large distances, intermediate range spin-dependent attraction 
dominated by two-pion exchange, and a short-range repulsion. One popular 
NN interaction that fits the experimental phase shifts well is Argonne 
At nuclear densities the two-body force is combined to a three-body force 
that is essential to reproduce the spectrum of light nuclei.'^ At very low 
densities, though, as found in neutron star crusts or the exterior of neutron- 
rich nuclei, the scattering length (a = —18.5 fm) and effective range (re = 
2.7 fm ) are most crucial to the physical properties of the system. The 
presence of a short-range repulsive core is important primarily in that it 
prevents collapse to a higher-density state. 

Schematically, Green's Function Monte Carlo projects out the lowest- 
energy eigenstate ^'o from a trial (variational) wave function "^v by treating 
the Schrodinger equation as a diffusion equation in imaginary time r and 
evolving the variational wave function up to large r. The ground state is 
arrived at using: 



evaluated as a branching random walk. The short-time propagator is 
sometimes taken from a Trotter-Suzuki approximation, but the exact two- 
nucleon propagator can also be employed. The ground-state energy £"0 can 





exp[-(i/ - Et)t]'1!v 
]Jexp[-(i/-i;T)AT]* 



(5) 
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be obtained from: 

^ {^v\H\^o) ^ {^o\H\^o) 
° (*y|vI/o) (vI/o|*o> • ^ ' 

In the case of more complicated interactions, one can use either the nu- 
clear Green's Function Monte Carlo methocP^or the Auxiliary Field Diffu- 
sion Monte Carlo approach!^ Schematically, the latter method reduces the 
spin-isospin dependence of the interaction operators from quadratic to lin- 
ear by means of the Hubbard-Stratonovich transformation. Given a generic 
operator O and a parameter A: 

V2ttJ 

As a result, the Monte Carlo sampling is no longer limited to the coordinate 
space positions of the particles, but now extends to sample auxiliary fields. 
Thus, the AFDMC algorithm limits the exponential growth of the spin- 
isospin states, and recovers the ground state in polynomial time. 

In these Quantum Monte Carlo superfluid simulations the trial wave 
function was taken to be of the Jastrow-BCS form with fixed particle num- 
berP 

* V = n /^(^«^) n /^(^■''^■') n fi^^MU '^^oi (») 

i¥'j id' i<j' 

and periodic boundary conditions. The primed (unprimed) indices corre- 
spond to spin-up (spin-down) neutrons. The pairing function 0(r) is a sum 
over the momenta compatible with the periodic boundary conditions. In 
the BCS theory the pairing function is: 

n " n 

The Jastrow part is usually taken from a lowest-order-constrained- 
variational methocP^l calculation described by a Schrodinger-like equation. 
The fixed-node approximation guarantees that the result for one set of pair- 
ing function parameters will be an upper bound to the true ground-state 
energy of the system. Variational results with the pairing function alone 
reproduce BCS calculations with finite particle-number projection. The 
parameters are optimized in the full QMC calculation, providing the best 
possible nodal surface, in the sense of lowest fixed-node energy, with the 
given form of trial function. Comparisons of fixed-node result^ in the cold 
atom system with recent AFQMC result^^^l indicate that the fixed node 
calculations are quite accurate. 
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In the case of closed-shell neutron drops the antisymmetric part of the 
trial wave function has the form 

i^{R,S)=\y^Det{4>^{n,s,)}\ , (10) 

where a = {n^j,mj} is the set of quantum numbers of single-particle or- 
bitals, and the summation of more determinants is done in order to have a 
trial wave function that is an eigenstate of and M. The single-particle 
basis is given by 

Mr,s) = $„,,(r) [l^.„Jf)6,™.(s)]^.„^ , (11) 

The radial components ^n.j are obtained solving the Hartree-Fock problem 
with a Skyrme force and their width is variationally optimized, Yi.mj are 
spherical harmonics and ^s.m^ are spinors in the usual up-down base. It 
is particularly important at low densities (low oscillator frequencies) to 
incorporate explicit BCS correlations in the trial wave function, essentially 
using single particle orbitals in the above equation in a form like Eq. [9j 

3. Infinite Matter Results 

The T — equations of state for homogeneous cold atoms and neutron 
matter of Ref. [8] are compared in Fig. [l] The horizontal axis is kpa, with 
the equivalent Fermi momentum kp for neutron matter shown along the 
top. The vertical axis is the ratio of the ground-state energy to the free 
Fermi gas energy (Efg) at the same density; as discussed in the previous 
section (Eq. Q), it must go to one at very low densities and decrease as 
the density increases and the interactions become important. The curve 
at lower densities shows the analytical result by Lee and Yang mentioned 
previously. The QMC results shown in this figure seem to agree with the 
trend implied by the Lee- Yang result. The neutron matter and cold atom 
equations of state are very similar even for densities where the effective 
range is comparable to the interparticle spacing. Near kpa = —10 the 
energy per particle is not too far from QMC calculations and measurements 
of the ratio ^ of the unitary gas energy to Epc shown as an arrow on the 
right (it corresponds to kpa = oo); previous calculations give ^ « 0.4, 
more recent results in Refs. 9|10 are slightly lower). At larger densities 



the cold-atom and neutron matter results start to diverge, due to: i) the 
neutron finite effective range, and ii) the fact that the neutron results also 
incorporate a simple attempt to include the 5 = 1, Ms — pairs. When 
the density is very low, the s-wave contribution is dominant so cold atoms 
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Fig. 1. Quantum Monte Carlo equation of state for cold atoms (squares) and neutron 
matter (circles). Also shown are the analytic expansion of the ground-state energy of a 
normal fluid (line) and an early Quantum Monte Carlo result at unitarity (arrow). 



and neutron matter agree very well. As the density increases, the effective 
range, as well as higher order terms in momentum and higher partial waves, 
become important. 

The pairing gap at T=0 is calculated using the odd-even staggering 
formula: 

A = [E{N + 1) - ^ [E{N) + E{N + 2)] ] (-1)^ , (12) 

where N is the number of particles. In Fig. [2] we plot the gap as a function 
of kpa for both cold atoms and neutron matter, taken from Ref. 8, BCS 
calculations are shown as solid lines, and QMC results are shown as points 
with error bars. QMC pairing gaps are shown from calculations of iV = 
66 — 68 particles. For very weak coupling, —kpa « 1, the pairing gap is 
expected to be reduced from the BCS value by the polarization corrections 
calculated by Gorkov and Mehk-Barkhudarov, A/Abcs — (l/4e)^/^, as 
mentioned in the previous section (Eq. ([s])). The QMC calculations at the 
lowest density, kpa = —1, are roughly consistent with this reduction from 
the BCS value. At slightly larger yet still small densities, where —kpa = 
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0(1) but kpVe « 1 for neutron matter, one would expect the pairing gap 
to be similar for cold atoms and neutron matter. The results at kpa = —2.5, 
where kpre ~ 0.35, support this expectation. Beyond that density the 
effective range becomes important and both the BCS and the QMC results 
are significantly reduced in relation to cold atoms where re ~ 0. 

In cold atoms, the suppression from BCS is reduced as the density 
increases, with a smoothly increasing fraction of the BCS results as we 
move from the BCS to the BEC regime. At unitarity the measured pairing 
gap^i^Hlil are 0.45(0.05) of the Fermi energy, for a ratio A/Abcs ~ 0.65, 
in agreement with predictions by QMC methodsEEE In the BEC regime 
where two fermions are tightly bound, the BCS and QMC values would 
both give a gap of half the binding energy of the pair. In neutron matter, 
the finite range of the potential reduces A/Ep as the density increases. We 
find a ratio A/Abcs that increases shghtly from |/ci?a| = 1 to 2.5, but then 
remains roughly constant up to \kpa\ — 10. 

In Refs. |36|41| new QMC values for neutron matter were compared to se- 
lected previous results. ^^ ^'^ ■^^ ■^^ The results of the QMC calculations 
are much larger than most diagrammatic approaches. As these approaches 



March 2, 2013 15:36 World Scientific Review Volume - 9in x 6in CGG Review 



10 J. Carlson, S. Gandolfi, and A. Gezerlis 

1.1 
1 

'|0.9 

mO.8 

0.7 
1.1 

1 

I 0.9 
UO.8 
0.7 

' 10 20 30 40 50 60 

N 

Fig. 3. Energy of neutrons confined in a harmonic well with frequency oj = 5 and 10 
MeV. Solid points are ab-initio calculations computed using AFDMC (red circles) and 
GFMC (blue squares). Open symbols and solid black lines represent results given by 
selected Skyrme forces. The violet solid lines show results obtained with the adjusted 
SLy4. 

assume a well-defined Fermi surface or calculate polarization corrections 
based on single-particle excitations it is not clear how well they can de- 
scribe neutron matter in the strongly paired regime, or the similar pairing 
found in cold atoms. Finally, the QMC results seem to qualitatively agree 
(at least for the lowest densities considered) with a determinantal Quantum 
Monte Carlo lattice calculation.!^ 

4. Neutron drops 

A neutron drop is an idealized system where neutrons, interacting via re- 
alistic nuclear forces, are confined in external potentials. They provide a 
simple model to study neutron-rich nuclei. In Refs. [42|43| the neutron-rich 
isotopes of oxygen and calcium have been studied as neutrons confined in 
external fields, and this model provided good results for both bound and 
excited states. In addition, these models can be used to constrain den- 
sity functionals in the large isospin-asymmetry region. Ab initio study of 
neutrons confined in different geometries showed that some Skyrme forces, 
typically fitted to nuclei with small isospin-asymmetry, are not accurate 
when dealing with pure neutron systems. The energies computed using 
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Fig. 4. Pairing gap of neutrons confined in a harmonic well with frequencies co = 5 (left) 
and 10 (right) MeV. Results were arrived at using different Hamiltonians with only a 
two-body force, and including different three-body forces are presented. 



QMC show that Skyrme forces tend to overbind confined neutrons, the 
main contribution to the difi^erence being the poorly constrained gradient 
term. In Fig. [3] we plot the energy of confined neutrons computed using 
QMC techniques, compared with selected Skyrme forces. 

The harmonic oscillator potential introduces another element into the 
pairing of neutrons, namely the mean field or the spacing of single-particle 
levels in the neutron drop. By studying the odd-even staggering of neutron 
drops with different frequencies, it is possible to examine the pairing as 
it evolves from small finite systems to infinite matter. The drops we have 
considered to date have harmonic oscillator frequencies of 5 and 10 MeV. In 
both these cases the pairing is much smaller than the single-particle spacing, 
so the BCS wave function is largely that occuring in the lowest shell of the 
drop. The GFMC calculations have included this pairing into the trial 
wave function, and have been used to examine both trap frequencies. At 
a frequency of 10 MeV, the pairing is reduced particularly for the larger 
drops. Early results are shown in Fig. [4] where the gap of confined neutrons 
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is calculated using different Hamiltonians. The dominant physics is 
provided by the NN interaction, with modest changes introduced by the 
three-nucleon interaction. Further studies of these systems are underway. 



5. Summary and Future Work 

In summary, this review has examined the application of Quantum Monte 
Carlo methods to strongly-paired matter, including homogeneous and 
trapped systems of atoms or neutrons. For infinite matter, the calculated 
equation of state and pairing gap match smoothly with the known analytic 
results at low densities, and provide important constraints in the strong- 
coupling regime at large kpa. The low-density equation of state can help 
constrain Skyrme mean-field models of finite nuclei. The pairing gap for 
low-density neutron matter is relevant to Skyrme-Hartree-Fock-Bogoliubov 
calculation^2ll of neutron-rich nuclei and to neutron-star physics, since it is 
expected to influence the behavior of the crust .'^ Similarly, results for neu- 
tron drops can also be used to constrain Skyrme and other energy-density 
functional approaches. Furthermore, the newly determined value of the gap 
implies that a new mechanism that makes use of superfluid phonons is com- 
petitive to the heat conduction by electrons in magnetized neutron stars 
Another consequence of the gap magnitude is related both to neutron-star 
observations and heavy-nuclei phenomenology: polarized neutron matter 
may be plausible within the context of magnetars, and has recently been 
attacked using Quantum Monte CarloE^E?! 

Microscopic many-body simulations will undoubtedly continue to strad- 
dle the divide between atomic and nuclear physics. Such simulations started 
with two species and equal populations, soon thereafter moving to spin- 
polarization. Further examination of the evolution of pairing from small to 
large systems is an intriguing area of study, including both cold atoms and 
neutron matter. For example, one avenue of future research is related to op- 
tical lattice experiments with cold atoms: to first approximation these are 
equivalent to periodic external potentials. In the nuclear case, an external 
potential would allow us to study the static response of neutron matter and 
would also facilitate the understanding of the impact on neutron pairing of 
the ion lattice that exists in a neutron star crust. Such microscopic results 
for the static response could provide further constraints on energy-density 
functionals used to describe the crust of neutron stars. 



March 2, 2013 15:36 World Scientific Review Volume - 9in x 6in CGG'Review 



Superfluid Pairing in Neutrons and Cold Atoms 13 

Acknowledgments This work was supported by DOE Grant Nos. DE- 
FG02-97ER41014 and DE-AC52-06NA25396. The computations shown 
were performed at the National Energy Research Scientific Computing Cen- 
ter (NERSC) and through Los Alamos Supercomputing. 



References 

1. S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 
(2008). 

2. B. Marcelis, B. Verhaar, and S. Kokkelmans, Phys. Rev. Lett. 100, 153201 
(2008). 

3. L. Luo and J. E. Thomas, J. Low Temp. Phys. 154, 1 (2009). 

4. N. Navon, S. Nascimbene, F. Chevy, and C. Salomon, Science 328, 729 

(2010) . 

5. J. Carlson, S.-Y. Chang, V. R. Pandharipande, and K. E. Schmidt, Phys. 
Rev. Lett. 91, 050401 (2003). 

6. G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, Phys. Rev. 
Lett. 93, 200404 (2004). 

7. J. Carlson and S. Reddy Phys. Rev. Lett. 95, 060401 (2005). 

8. A. Gezerlis and J. Carlson, Phys. Rev. C 77, 032801(R) (2008). 

9. M. M. Forbes, S. Gandolfi, and A. Gezerlis, Phys. Rev. Lett. 106, 235303 

(2011) . 

10. J. Carlson, S. Gandolfi, K. E. Schmidt, and S. Zhang, Phys. Rev. A 84, 
061602(R) (2011). 

11. M. J. H. Ku, A. T. Sommer, L. W. Cheuk, and M. W. Zwierlein, Science 
335, 563 (2012). 

12. Y. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, Nature 451, 689 
(2008). 

13. J. Carlson and S. Reddy, Phys. Rev. Lett. 100, 150403 (2008). 

14. A. Schirotzek, Y. Shin, C. H. Schunck, and W. Ketterle, Phys. Rev. Lett. 
101, 140403 (2008). 

15. U. Lombardo and H.-J. Schulze, Lecture Notes in Physics ( Springer- Verlag, 
Berlin, 2001), Vol. 578, p. 30. 

16. D. J. Dean and M. Hjorth- Jensen, Rev. Mod. Phys. 75, 607 (2003). 

17. S. Gandolfi, J. Carlson and Steven C. Pieper, Phys. Rev. Lett. 106, 012501 
(2011). 

18. D. Page, J. M. Lattimer, M. Prakash, and A. W. Steiner, Astrophys. J. 707, 
1131 (2009). 



19. D. Page, S. Reddy |arXiv:1201.5602 (2012) 



20. D. N. Aguilera, V. Cirigliano, J. A. Pons, S. Reddy, and R. Sharma, Phys. 
Rev. Lett. 102, 091101 (2009). 

21. V. Cirigliano, S. Reddy, and R. Sharma, Phys. Rev. C 84, 045809 (2011). 

22. N. Chamel, S. Goriely and J.M. Pearson, Nucl. Phys. A812, 72 (2008). 

23. T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957). 



March 2, 2013 15:36 World Scientific Review Volume - 9in x 6in CGG Review 



14 J. Carlson, S. Gandolfi, and A. Gezerlis 

24. L. P. Gorkov and T. K. Mclik-Barkhudarov, JETP, 40, 1452 (1961) [Soviet 
Phys. JETP 13, 1018 (1961)]. 

25. H. -J. Schulze, A. Polls, and A. Ramos, Phys. Rev. C 63, 044310 (2001). 

26. J. M. C. Chen, J. W. Clark, R. D. Dave, and V. V. Khodel, Nucl. Phys. 
A555, 59 (1993). 

27. J. Wambach, T. L. Ainsworth, and D. Pines, Nucl. Phys. A555, 128 (1993). 

28. H.-J. Schulze, J. Cugnon, A. Lojeune, M. Baldo, and U. Lombardo, Phys. 
Lett. B375, 1 (1996). 

29. A. Schwenk, B. Friman, and G. E. Brown, Nucl. Phys. A713, 191 (2003). 

30. H. Miither and W. H. Dickhoff, Phys. Rev. C 72, 054313 (2005). 

31. A. Fabrocini, S. Fantoni, A. Y. Illarionov, and K. E. Schmidt, Phys. Rev. 
Lett. 95, 192501 (2005). 

32. L. G. Cao. U. Lombardo, and P. Schuck, Phys. Rev. C 74, 064301 (2006). 

33. J. Margueron, H. Sagawa, and K. Hagino, Phys. Rev. C 77, 054309 (2008). 

34. S. Gandolfi, A. Yu. Illarionov, S. Fantoni, F. Pederiva, and K. E. Schmidt, 
Phys. Rev. Lett. 101, 132501 (2008). 

35. T. Abe and R. Scki, Phys. Rev. C 79, 054002 (2009). 

36. S. Gandolfi, A. Yu. Illarionov, F. Pederiva, K. E. Schmidt, and S. Fantoni, 
Phys. Rev. C 80, 045802 (2009). 

37. R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C 51, 38 (1995). 

38. Steven C. Piepor, V. R. Pandharipande, R. B. Wiringa and J. Carlson, Phys. 
Rev. C 64, 014001 (2001). 

39. K. E. Schmidt and S. Fantoni, Phys. Lett. B446, 99 (1999). 

40. V. R. Pandharipande and H. A. Bethe, Phys. Rev. C 7, 1312 (1973). 

41. A. Gezerlis and J. Carlson, Phys. Rev. C 81, 025803 (2010). 

42. S. Gandolfi, F. Pederiva, S. Fantoni and K. E. Schmidt, Phys. Rev. C 73, 
044304 (2006). 

43. S. Gandolfi, F. Pederiva and S. a Boccara, Eur. Phys. J. A 35, 207 (2008). 

44. J. Carlson, S. Gandolfi, Pieter Maris, James Vary and Steven C. Pieper, in 
preparation. 

45. A. Gezerlis, Phys. Rev. C 83, 065801 (2011). 

46. A. Gezerlis and R. Sharma, Phys. Rev. C 85, 015806 (2012). 



